Pygmy Quadrupole Resonance in Skin Nuclei 
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The electric quadrupole response is investigated theoretically by HFB and QPM calculations along 
the Sn isotopic chain with special emphasis on excitations above the first collective state and below 
the particle threshold. Depending on the asymmetry, additional strength clustering as a group of 
states similar to the known PDR mode is found. The spectral distributions and electric response 
functions are discussed. The transition densities of these Pygmy Quadrupole Resonances (PQR) 
states are closely related to the neutron excess and showing special features being compatible with 
an oscillation of the neutron skin against the nuclear core. These features may indicate a new 
nuclear quadrupole mode connected to the skin configuration. 



INTRODUCTION 



The new radioactive beam facilities allow to probe 
hitherto unaccessible regions of nuclear masses far from 
the valley of stability The prevalence of one type of 
nucleons in medium and heavy nuclei may result in the 
formation of a nuclear skin, containing either pure pro- 
ton or neutron matter, respectively. After the detection 
of halos in light nuclei 0, Q , the skin phenomena adds 
another new aspect to the dynamics of strongly asym- 
metric nuclear matter, observable in medium and heavy 
exotic nuclei. Besides the genuine importance for nu- 
clear structure physics, skin nuclei are of particular rel- 
evance for nuclear astrophysics 0, H[. For example, the 
existence of low-energy excitation modes related to the 
neutron skins could influence the rapid (r) neutron cap- 
ture process in low-energy (n,7) reactions and, therefore, 
will affect directly the synthesis of heavy nuclei in stellar 
environments. 

One of the most interesting results was the discovery 
of a new dipole mode at energies below and close to the 
particle emission threshold (y, observed in stable and 
unstable nuclei with isospin asymmetry N/Z > 1. Typ- 
ically that strength is found as a bunching of discrete 
l - states with very similar spectroscopic features. Al- 
though located below the particle threshold, the clus- 
tering of states in a narrow energy region resembles a 
resonance phenomena. This extra dipole strength was 
unexpected and could not be explained as just belong- 
ing to the low-energy part of the giant dipole resonance 
(GDR) Q. Because of the unique character of the ex- 
citation mechanism the mode was named pygmy dipole 
resonance (PDR). The PDR is of high current interest 
and its properties are investigated by a larg e number 
of theoretical and experimental groups |6l-[l9|. Most of 
the approaches agree on the hypothesis that the PDR 
is an isospin-mixed mode due to excitations of identi- 
cal particles, neutrons (protons) from the nuclear sur- 
face, corresponding to an oscillatory motion of the nu- 
clear skin against the core [13|-|19(. Despite remaining 
differences among the theoretical approaches (in particu- 
lar between relativistic and non-relativistic calculations), 



there is agreement on the gross features of these excita- 
tions as a new nuclear mode, being closely related to and 
correlated with the size of the nuclear skin, as discussed 
in detail e.g. in refs. [H, [H, [H . 

An obvious question, coming up immediately in this 
connection, is to what extent the presence of a neu- 
tron (proton) skin will affect also excitations of other 
multipolarities. Promising candidates are low-energy 2 + 
states, in excess of the spectral distributions known from 
stable nuclei close to the valley of stability. Previous 
theoretical studies are found in [21M23I ] . Similar results 
were presented in our recent quasiparticle random phase 
approximation (QRPA) calculations in 120 Sn 24] where 
low-energy electric quadrupole strength, located well be- 
low the isoscalar giant quadrupole resonance (ISGQR) 
25] was observed. In order to clarify whether this is 
a general feature of the skin nuclei we have studies the 
quadrupole response along the Sn-isotopic chain. This re- 
port contains results based on a Hartree-Fock-Bogoljubov 
(HFB) description of the ground state (2(| by usinga phe - 
nomenological energy-density functional (EDF) 0, [l3j . 
The excited states are calculated with QRPA theory, 
using the multiphonon description of the Quasiparticle- 
Phonon model (QPM) (f^. 



II. THEORETICAL MODEL 

The building blocks of the model are given by the 
quasiparticle states obtained from a HFB description of 
the nuclear ground states and coherent superpositions of 
two-quasiparticle states treated by QRPA theory. A de- 
tailed description of the approach is found in 0, [EH , here 
we give only a brief overview. The approach is based on 
the Hamiltonian 

H = Hmf + H res , (2.1) 

The mean- field part Hmf defines the single particle prop- 
erties and, as such, accounts for the ground state dynam- 
ics, including potentials and pairing interactions for pro- 
tons and neutrons, respectively. That part is related to 
an EDF 0, HH constructed such that also dynamical ef- 
fects beyond mean-field can be taken into account. That 
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FIG. 1. EDF and HFB calculations of nuclear skin thicknesses 
Sr in comparison with skin thicknesses derived from charge 
exchange reactions by Krasznahorkay et al. 29,30]. The skin 
thickness is defined with the deference between neutron and 
proton rms radii with the equation Sr — \/ < r% > — v/< r| >. 



goal is achieved in practice by using fully microscopic 
HFB potentials and pairing fields as input but performing 
a second step variation with scaled auxiliary potentials 
and pairing fields readjusted in a self-consistent manner 
such that nuclear binding energies and other ground state 
properties of relevance are closely reproduced. 

The residual interactions contained in H res are treated 
separately. After having derived the ground state single 
quasiparticle spectra and wave functions the nuclear ex- 
cited states are described by QRPA phonons of various 
multipolarities. Following the notation of )27j, the state 
operators are defined as a linear combination of two- 
quasiparticle creation and annihilation operators corre- 
sponding to excitations propagating forward and back- 
ward in time, respectively: 
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(2.2) 



Here, j = (nljmr) is a single-particle proton or neu- 
tron state; At and are proton and neutron two- 



quasiparticle creation and annihilation operators, respec- 
tively, with configuration amplitudes given by the time- 
forward and time-backward amplitudes i/jjj and (p^ l i 
27] , respectively. The operators are coupled to total 
angular momentum A with projection \i by means of 
the Clebsch-Gordan coefficients C-^., m , = (jmj'm' |A/z). 
The phonons obey the QRPA equation of motion 



(2.3) 



accounting now also for the residual interactions in the 
particle-hole and the particle-particle channels, respec- 
tively. By solving the eigenvalue problem, eq. 12. 31 we 
obtain the phonon excitation energies E a and the state 
vectors, eq. 12.21 

In more detail, the residual interactions are given by 



TT _ TfPh 

n r€.s — M 



H 



ph 
SM 



ttVV 

M 



(2.4) 



At present, we do not attempt a fully self-consistent 
derivation of H res but follow the empirical approach of 



1* 



— * — exp 
— ■— QRPA 



£- 0,0 I ' ' ' ' ' ' ' ' ' ■ ' ' ' ' ' 

m 104 108 112 116 120 124 128 132 



FIG. 2. QRPA calculations of the energies (upper panel) 
and transition strengths of the first 2^ states in 104 ~ 134 Sn 
compared to experimental data [3214391 ]. 
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FIG. 3. Calculations the neutron and proton contents in the 
structure of the 2j" states in terms of weight factors ujj 1 j 2 (A/xi 
where the sum is taken over the neutron or proton components 
entering in the structure of a phonon. 



the QPM. Hence, we use separable multipole-multipolc 



H M and spin-multipole Hg M interactions both of 
isoscalar and isovector type in the particle-hole (ph) and 
multipole pairing iJfJ in the particle-particle (pp) chan- 
nels, respectively [3lj . Further details for the definition 
of the interactions and the determination of the coupling 
str eng th parameters have been discussed extensively e.g. 
in [13j and it will not be repeated here. We only mention 
the interesting observation that in the nuclei considered 
here the isoscalar quadrupole constant k is rather inde- 
pendent of the nuclear mass with only slight variations, 
typically in the order of 1% between neighboring stable 
nuclei and less than 5% outside the valley of stability. 
This allows safe extrapolations into the hitherto unex- 
plored regions toward 100 Sn and beyond 132 Sn. 
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III. QUADRUPOLE EXCITATIONS IN SN 
ISOTOPES 

A. The quadrupole response functions 

Before turning to the discussion of the quadrupole re- 
sponse, it is worthwhile to refer again to our previous 
studies. The HFB calculations [l3[ have shown that in 
the mass region A>106, the tin nuclei exhibit a neutron 
skin while the situation reverses for the lower masses to 
a tiny proton skin in 100-104 Sn, as also seen in Fig. 1. 
A neutron skin, being a static property of the nuclear 
ground state, manifests itself also in the dynamic sector 
in terms of the low-energy PDR dipole excitations. Ac- 
cording to theory, the phenomenon is quite ubiquitous, 
shown e.g. by our previous studies in the tin isotopes 
and the N=50, 82 nuclei @,[ll|. Meanwhile, also 
a growing evidence from experiments is obtained. In gen- 
eral, the PDR states are of mixed isospin symmetry and 
their total transition strength is closely related to the size 
of the neutron skin [l3j . 

Now coming to the present investigations, we first dis- 
cuss the properties of the traditionally known parts of the 
quadrupole spectrum. These considerations serve also as 
tests before focusing on a new sector of the quadrupole 
spectrum. The QRPA results for energies and B(E2) 
transition probabilities of the 2^ state are displayed and 
compared to data [34j - t36| in Fig. 2. The excitation ener- 
gies are remarkably stable over the mass range, except at 
the shell closure at A=132. Also the transition strengths 
vary only mildly in a narrow band for most of the isotopes 
but decreasing when 132 Sn is approached. In 104 Sn we 
predict a 2f state with excitation energy B x =1.221 MeV 
and B(E2;g.s.-> 2+)=0.10 e 2 b 2 . Thus, we confirm the 
experimentally observed tendency of decreasing B(E2) 
values (also found in shell model calculations) in nuclei 
with extreme proton-to- neutron ratios N/Z toward 100 Sn 


Here, the region around the double magic nucleus 
132 Sn is of particular interest. The restoration of 
N=82 neutron-shell gap was recently confirmed by high- 
precision mass measurements (40j . From a systematic 
study of 2^ states the authors of [4l| found that the 
product of transition strength and excitation energy is 
almost constant, i.e. B(E2; g.s.->2f )w 1/E(2^). First of 
all, this means that the contributions of the 2+ states to 
the electromagnetic energy weighted sum rule (EWSR) is 
to a good approximation a constant. Since the EWSR is 
a conserved quantity, this puts constraints on the higher 
lying energy weighted strength which, consequently, is 
separately conserved as well. 

For the explanation of this effect the collectivity of the 
excited states plays a crucial role. Microscopically, collec- 
tivity is directly related to the coherent, i.e. in-phase, ex- 
citation of a large number of two-quasiparticle pairs with 
respect to a transition operator. Favorable candidates are 
those states which are shifted by interactions into spec- 
tral regions away from the poles of the non-interacting 



two-quasiparticle Green function (i.e. the 4-point func- 
tion). Namely, more collective states have lower excita- 
tion energies and respectively stronger transition rates 
(because the nuclcons contribute coherently). An excep- 
tion from this rule are nuclei around the shell closures as 
observed in N=126 (Pb region) [38j and N=82 41]. Our 
QRPA results for 130_134 Sn agree with the conclusions 
drawn in ref . 0, HU . 

The abnormal behavior of the B(E2) in nuclei around 
132 Sn (see Fig. 2) can be traced back to an increase of 
the amplitude of the lgg/22d 5 / 2 two-quasiparticle proton 
component in the first QRPA 2+ state in 132 Sn, amount- 
ing for about 17% of the state vector. Compared to the 
neighboring nuclei 130 Sn and 134 Sn, where this compo- 
nent contributes with probabilities of 1.5% and 0.9%, re- 
spectively, this is a substantial reordering of strength. 
Results of the summed proton/neutron amplitudes in 
terms of weight factors 

= E \& - i^i 2 * 10 °[ % ]' 

contributing to the state vector of the 2+ state in Sn 
isotopes are presented in Fig. 3. 

After having fixed accurately the model parameters 
and since both energies and B(E2)-values of the lowest 
quadrupole states are well reproduced, we expect that 
the calculation will lead to reliable results also for the 
quadrupole states at intermediate excitation energies. 
Our particular interest is to investigate whether there 
is a correlation between the neutron (proton) skin and 
quadrupole states in the energy region below and around 
the particle separation threshold, similar to the dipole 
case. For completeness and as further test cases we also 
include the higher-lying collective states up to the isovec- 
tor giant quadrupole resonances (IVGQR) (2f| and above 
until E x = 35 MeV. Similarly to the QRPA results for 
the PDR modes, in Sn isotopes (N= 66-84) we observe 
a sequence of quadrupole states with common spectro- 
scopic features, located just above the collective 2^ and 
below the ISGQR. The lowest-lying of them, located typ- 
ically in the energy region E x w 2-4 MeV (with exception 
the lighter tins and 132 Sn), are non-collective, with pre- 
dominantly neutron structure (the proton contribution is 
less than a few percents). With increase of the excitation 
energy above E x s» 4 MeV toward the particle emission 
threshold, the collectivity of the states increases, though 
they have still more pronounced neutron components for 
A>106. 

We suggest that the low-energy quadrupole excita- 
tions with predominantly neutron structure in tin iso- 
topes could belong to a new quadrupole mode, namely, 
a neutron PQR which manifests itself as a neutron sur- 
face vibrations. A supporting argument in this direction 
could be a possible connection between the B(E2) transi- 
tion rates of the PQR states and the neutron number of 
the nucleus related to the thickness of the neutron skin 

From systematic studies of B(E2) transitions of 
quadrupole states with excitation energies E x sa 2-4 MeV 
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FIG. 4. Calculations of isoscalar (a) and isovector (b) electric quadrupole strengths in comparison with the total QRPA B(E2) 
transition strengths (c) up to 35 MeV in Sn isotopes. 
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FIG. 5. QRPA proton (dash line) and neutron (solid line) quadrupole transition densities summed over 2 + excited states in a 
given energy region (indicated over every plot) in 104 - 120 > 134 s n isotopes. 
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in the tin isotopes we find a correlation between the 
summed over non-collective neutron modes B(E2) tran- 
sitions and the number of the neutrons of the nucleus 
related to the neutron skin thickness. Theoretical cal- 
culations of B(E2) spectral transition strength distribu- 
tions in 104 . 12 o. 13 4g n isotopes are presented in Fig. |4] 
A sizable effect of increased B(E2) at E x « 2-4 MeV is 
observed for the heaviest tin isotopes - 130 Sn and 134 Sn. 
For 132 Sn the neutron excitations in this energy range are 
suppressed due to the N=82 shell closure. In the lighter 
tin isotopes the reduced neutron number decreases the 
collectivity of the 2+ state which becomes non-collective 
in 104 Sn. At the same time the proton contribution to 
the B(E2) strength located in the energy range 2-4 MeV 
increases toward 104 Sn and brings to more intensive pro- 
ton quadrupole excitations in 104 Sn there (see Fig. 
Consequently, the 104 Sn nucleus appears to be an op- 
posing point where the neutron skin chances to proton 
and the neutron PQR to proton PQR, respectively. The 
investigations of the nature of the PQR excitations in 
104- 134 g n nuc i e j r eveal that these states are predomi- 
nantly of isoscalar character in the heavier tin isotopes. 
The contribution of isovector strength in 130 Sn and 134 Sn 
in the PQR region is negligibly small but it increases to- 
ward lighter tin isotopes due to the larger proton content 
of the state vectors. The calculations of the distribu- 
tions of the isoscalar Mo and isovector Mi quadrupole 
strengths in 104 ~ 134 Sn up 35 MeV are presented in Fig. 
Qb>,b, respectively. The isoscalar (isovector) electric ma- 
trix elements Mgm(A^) for A = 2 are defined with the 
equation: 



M 0(1) (2+) 



p 

E 



rgy a „(n fc ) 



± £- 2 

k 



g.s. 



(3-2) 

The sign '+' in this equation corresponds to the 
isoscalar case for which the entering 2 + states have the 
same signs of the proton and neutron summed single- 
particle matrix elements. The sign '-' is related to the 
isovector case where the proton and neutron summed 
single-particle matrix elements are with opposite signs, 
respectively. 



B. Quadrupole Transition Densities 

The characteristics of the PQR states are seen most 
clearly by inspection of transition densities, defining the 
non-diagonal parts of the one-body density matrix [l3| . 
In this way, we expect a deeper insight into the spatial 
pattern and location of the transition strength, similar 
to the analysis for dipole excitations [l3j]. In Fig. [5] 
results for the proton and neutron transition densities 
of the isoscalar 2+ state in several tin isotopes are dis- 
played. From the QRPA state vectors one finds that the 
2f state changes from collective in the middle of the shell 
(between N=50 and N=82) to non-collective toward the 
shell closure. 



In 106_134 Sn most of the quadrupole states below the 
particle emission threshold are of isoscalar character. 
However, some mixed symmetry configurations, indi- 
cated by the opposite signs of the proton and neutron 
amplitudes, are also observed in this region, as seen in 
the third 120 Sn panel, FigJS] At lower energies, E x =2-4 
MeV and E x =0. 8-3.7 MeV in 120 Sn and 134 Sn, respec- 
tively the main part of the oscillations is coming from 
less strongly bound neutrons, forming a skin-like surface 
layer [l3j], located predominantly at the nuclear periph- 
ery, where the nuclear radii are extended above an aver- 
age value of R~5 fm up to RwlO fm in 120 Sn and up to 
R«20 in 134 Sn. At the same time we find a very small 
proton contribution in this space region. 

For the double-magic 132 Sn the first QRPA 2+ (E x =3.9 
MeV) and the second QRPA 2+ (E x = 5.3 MeV) states 
exhibit properties similar to the neutron states discussed 
in the other tin isotopes in the same energy range. These 
states incorporate strength related to oscillations of neu- 
trons at the surface region. In addition in the second 
2 + state in 132 Sn in the nuclear interior (R < 5 fm) the 
protons and neutrons vibrate in-phase where proton an 
neutron parts are of comparable amplitude. 

With change in the structure of the quadrupole states 
at E x —2- 5 MeV in respect of increase in the proton 
contribution in the states vectors toward 104 Sn the pro- 
ton and neutron transition densities undergo a change in 
their behavior as it is demonstrated in Fig. [SJ The pro- 
ton oscillations at the nuclear surface at E x = 2- 4.3 MeV 
in Sn prevail over the neutron oscillations. The ob- 
servations of the change of the behavior of the transition 
densities in 104 Sn could be a sign for transformation of 
the character of the PQR from neutron to proton (simi- 
larly to the PDR case jll [TBI fl9j). 

With the increase of the excitation energy E x toward 
the particle separation threshold, the quadrupole states 
become more collective with a larger admixture of isovec- 
tor components to the state vectors. Nevertheless, some 
isospin effects are still present up to the particle emis- 
sion threshold (even strongly hindered). In this case, 
in 106 - 130 Sn isotopes we observe more isovector type of 
oscillations inside the nucleus while at the surface neu- 
trons dominate. The nature of these higher-lying states 
to our opinion could not be related directly to skin ef- 
fects as they collect a lot of strength (some of the them 
have transition probability B(E2)t of order of 10 W.u., 
which is of the same size as for the collective 2+ and the 
IS(V)GQRs) and the contribution of nucleons from the 
nuclear interior is significant. At the same time the role 
of the isospin effects, in particular the neutron skin for 
energies close to the threshold has to be further inves- 
tigated. The last two plots (down) in every tin nucleus 
presented in Fig. [S] show the proton and neutron transi- 
tion densities of the ISGQR and IVGQR, respectively. 

Overall, the transition densities show an interesting 
dependence on the excitation energy: in all nuclei the 
isoscalar character of the lowest 2 + state and the GQR 
is evident from the in-phase behavior of proton and neu- 
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tron components, while proton and neutron transition 
densities carry opposite phases in the isovector GQR re- 
gion. Together with the almost identical radial shapes 
these results follow almost perfectly well the generally 
accepted rules for collective quadrupole states. They are 
well interpreted as surface and volume oscillations of a 
liquid drop [43j | . However, the transition densities re- 
lated to the additional quadrupole states which may be 
considered as PQR excitations are not following any of 
the known rules but are showing quite unusual proper- 
ties. As seen from Fig. [5]the proton components continue 
to behave like vibrational variations of the nuclear den- 
sity radius while the neutron transition densities develop 
a rather different nodal pattern. The nodal structures 
correspond to processes where a (tiny) portion of nuclear 
matter is shuffled around the nuclear radius, leaving the 
latter almost unaffected, as indicated by the radial node 
occurring at or close to the nuclear half density radius, 
R ~ 5 fm. In a liquid drop picture, these features can 
be interpreted as vibrational excitations of the nuclear 
surface given by a superposition of variations in position 
and diffuseness of the nuclear surface. The latter mode 
was already anticipated by Bohr and Mottelsen al- 
though hitherto experimentally never verified, at least 
not in stable nuclei. Our results give evidence that such 
modes may indeed exist in charge-asymmetric nuclei. 

This interpretation is supported by a result of Pethick 
and Ravcnhall who derived in ref. [45j a relation con- 
necting the thickness of the nuclear skin to the sur- 
face tension. According to their result, the surface ten- 
sion decreases with increasing neutron excess. This is a 
statement concerning the static properties of the nuclear 
ground state density which is at least qualitatively fully 
supported by our microscopic HFB ground state, densi- 
ties. Now, as mentioned, the shapes of the PQR neutron 
transition densities are compatible with a reordering of 
neutron skin matter around the nuclear surface, hence 
inducing a local variation of the neutron excess in the 
surface region, resulting in oscillatory alterations of the 



surface tension. 



IV. CONCLUSIONS 

Quadrupole states were investigated in 104 ~ 134 Sn nu- 
clei up to excitation energies of 35 MeV. From QRPA 
calculations of states, the QPM model parameters are 
reliable determined. For the double- magic 132 Sn an in- 
crease of the B(E2; g.s. — > 2f ) compared to the neighbor 
130 Sn and 134 Sn is observed and explained. For 104 Sn, the 
first 2 + state is predicted and its spectroscopic proper- 
ties are determined. From the analysis of isoscalar and 
isovector electric quadrupole strength distributions the 
nature of the quadrupole states is investigated and a sep- 
aration between collective isoscalar, isovector and low- 
energy mixed symmetry states is achieved. The spectral 
distribution of the low-energy mixed symmetry states in 
io4-i34g n are c i us tering in a confined energy region such 
that they may be considered forming a PQR. An increase 
of the lowest-energy non-collective strength toward 134 Sn 
due to mainly neutron excitations is found. In addi- 
tion, the correlation of the PQR transition strength with 
the neutron (proton) skin thickness manifests itself via 
a transition from a neutron PQR to a proton PQR in 
104 Sn, the mass region where the neutron skin reverses 
into a proton skin. By means of quadrupole transition 
densities the nature of the PQR is clarified. In general, 
the PQR resembles the properties of the PDR and could 
be related to skin oscillations of one type of nucleons. 
Even though, with the increase of the collectivity of the 
quadrupole excitations the reliable distinction of the pure 
skin neutron (proton) quadrupole oscillations seems to be 
a difficult task, especially for nuclei with small or mod- 
erate neutron excess. In this respect, nuclei with larger 
isospin asymmetry like 130 Sn and 134 Sn could be a good 
candidates for experimentally explorations. 
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